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Abstract 

For thermal equilibrium systems it is shown, how the Kubo-Martin-Schwinger 
boundary condition may be used to factorize the generating functional of 
Green functions at least on the level of the full two-point function. Genuine 
non-equilibrium system exhibit correlations that one may also incorporate in 
the path integral. One one hand this provides a natural tool for a perturbative 
expansion including these correlations. On the other hand it allows to prove 
that in general non-equilibrium systems the generating functional does not 
factorize. 
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I. INTRODUCTION 

To a great extent, modern physics has turned its attention to non-equilibrium systems. 
Relativistic heavy-ion collisions JlJ and the early universe [|] are two very complementary 
examples for such systems. They have in common however, that their theoretical description 
to some extent follows guidelines of relativistic quantum field theory. 

Among these guidelines is the feature that causal propagation of particles requires empty 
(particle) states to propagate forward in time, whereas occupied states propagate backward 
in time. For the vacuum this guideline leads to the Feynman boundary condition of the 
two-point Green function (propagator), for thermal equilibrium states it leads to the Kubo- 
Martin-Schwinger (KMS) boundary condition [|3|. For non-equilibrium states however, this 
guideline requires to double the Hilbert space for the theoretical description: Since the 
occupation of states may change with time, one has to carry time-forward as well as time- 
backward boundary condition along throughout a computation. 
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Causality therefore imposes an additional 2x2 matrix structure on field theoretical 
descriptions of non-equilibrium states. This holds also in the path integral formulation of 
quantum field theory, which we consider of particular importance for some physical prob- 
lems: It is the basis of correlation expansion schemes [|j] as well as of derivations of finite 
temperature diagram rules Furthermore, it is used in lattice gauge theory, and may be 
found as the principal object in the recent discussion of the equivalence of different thermal 
equilibrium limits of quantum field theory J7LH . 

With the present paper we wish to demonstrate, that one may greatly reduce the calcula- 
tional effort required to take the Hilbert space doubling into account in path integrals. This 
is possible due to a diagonalization scheme for the 2x2 matrix structure . Apart from 
this more technical improvement (which might be considered conceptual as well), we also 
address the question of boundary conditions in general. We present a recipe, how to intro- 
duce initial correlations (which are a genuine non-equilibrium effect) into the path integral 
formalism [|TT|,|T^]. 

The doubling of the Hilbert space exists in several different flavors, the two most common 
are the Schwinger-Keldysh or closed time path method (CTP) [|13j and thermo field dynamics 
(TFD) ||14|| . Since one may show, that apart from conceptual differences they lead to identical 
results, we will use them as equivalent in the present paper: Computations are carried out 
in the CTP formalism, but we will make use of results obtained in TFD. 

The paper is organized as follows: In the next section, we introduce our technical im- 
provement into the path integral representation of a free quantum field. In section 3, the 
considerations are extended to the interacting case. Section 4 is devoted to the problem of 
initial correlations in terms of Wightman functions. In section 5 we replace the Wightman 
functions of the correlations by a cumulant expansion, and close the argumentation of this 
paper by relating this expansion to a classical integral equation including temporal boundary 
conditions. 



II. BOUNDARY CONDITIONS IN THE PATH INTEGRAL 

The first step in our discussion is to introduce the concept of diagonalization in the 
Schwinger-Keldysh or TFD method. To this end we study the physically trivial example of 
a free scalar quantum field. Throughout the paper we will use (in a somewhat sloppy fashion) 
two different notations for space-time arguments of fields, <p x = <j)(x). The Lagrangian for 
the free scalar field is 

C (<f>) = --<f> x (□ + m 2 ) (t> x . (1) 

The action integral extends over the whole three-dimensional coordinate space, but in time 
direction has a more complicated structure. Due to the fact, that our system has a pre- 
sumably irreversible time evolution, the time component of the action integral follows the 
contour depicted in fig. 1: Forward in time as well as backward in time, separated by an 
infinitesimal amount above and below the real time axis (For brevity, we do not discuss the 
foundation for this action integral, but rather refer to a review article on the doubling of 



the Hilbert space [[R|]). In terms of the time contour, the generating functional for Green 
functions is 
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and Green functions are contour ordered along C. Thus, in taking functional derivatives 
with respect to J it will make a difference whether the space-time arguments of J are on 
the upper or lower part of the contour C. 

Since under the contour ordering symbol the fields on the upper and lower branch of the 
contour commute at any time, we may treat them as two completely different fields. This is 
also mathematically more correct, because in thermal states the KMS condition induces a 
doubling of the Hilbert space [HJ. Indeed, with this step one recovers the concept of thermo 
field dynamics |[TJ]. The two fields will be distinguished by a lower index, <pi^ x is the field 
on the upper branch of the time contour. 

The generating functional for the Green functions of these two free fields is obtained as 



X>[01,0 2 ] exp 



i <Tx £ 



-om, 



+ <j>l,xJ\,x — 4>2,xJ2,. 



with the free Lagrangian 



£o(0i, 02) = -^0i,x ( D + rn 2 ) <f> ltX + -02,2 ( D + rn 2 ) 2 , 



(3) 



(4) 



and two independent classical source currents J\ and J 2 - 

Since we are dealing with free fields, the generating functional may be obtained in closed 
form, 



Zq [Ji, J2 



exp 



-- jd xd y[ J ljX Dj(x - y)Ji, y - J2,x.Dj(x - y)J 2 , y 



(5) 



The question now arises, which boundary conditions to chose for the two propagators Dj(x- 
y) and Dj(x — y). Both are resolvents to the free Klein-Gordon equation, e.g. 



(D x + m 2 ) Dj(x -y) = -5\x - y) 



(6) 



We may use Occam's razor as a guideline at this point: The simplest choice clearly are 
retarded and advanced boundary conditions, since they are free of any influence of the 
occupation of states. Due to the introduction of the path integral as running along a 
contour forward and backward in time, we therefore chose Dj(x — y) to be the free retarded 
propagator 



Dj(x -y) = D*(x -y) = -i &(x° - y°) [0(x), <j>{y)\ 



(7) 



The other propagator concurrently is chosen to have an advanced boundary condition in 
time, 



Dj{x -y) = D$(x -y) = -i Q(y° - x°) [<f>(x), <j>(y)] 



(8) 
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With this ansatz we have completely decoupled the two fields on the upper and lower time 
branch of the time-path depicted in fig. 1. In particular it is obvious, that mixed functional 
derivatives of this generating functional are zero when the sources are set to zero: 

Z [Ji,J2}) =0. (9) 

/ J 1 =J 2 =0 



5J 1 (x)5J 2 (y) 



One may now introduce an interaction, and expand the theory perturbatively in terms of 
the propagators Dj and Dj. The perturbation expansion is incredibly simple, since it is free 
of any occupation number factors, like e.g., Bose-Einstein distribution functions. Moreover, 
several types of diagrams in such a perturbative expansion are trivially zero, like e.g. a loop 
integral over two retarded propagators. 

However, it is also obvious that the sources J\ and J 2 introduced in the path integral are 
not the causal sources necessary to calculate physically meaningful quantities. Instead we 
need to define physical external sources j\ an d j 2 , which are obtained by a linear functional 
from the diagonal sources J\ and J 2 . The linearity assures, that the condition J± = J 2 = 
simply translates into j\ = j 2 = 0. 

We obtain them by demanding that the mixed second functional derivatives of the gen- 
erating functional fulfill the Kubo-Martin-Schwinger (KMS) boundary condition || 



Sji(x)Sj 2 (y) ^ l, ^' :j2 '/ 

Z \jl,j2]) ■ (10) 



Sj 2 (x° - i(3,x)8j 1 (y) 



31=32=0 



This is equivalent to a periodicity condition of the fields in the imaginary time direction. 

To obtain the linear functional relating the j-sources to the J-sources, it is convenient 
to switch to a Fourier representation 



oo 



= Iwf^ {e- iEX °^\E,k)+e^°t\E,k)) (11) 

o 

for % = 1,2, and similarly for the sources Jj. The shift in the imaginary time direction 
then simply translates into a multiplication with a Boltzmann factor. One possible way 
to express the KMS condition amounts to a linear relationship between the two types of 
sources, obtained from the Boltzmann factor as 

r 3 B(n B (E))r 3 , (12) 




jt\E,k) 
{ -\E,k)) \4-\E,k) 



T 



with a matrix 



B{n B W)=( 1 _i nB{E!) - nB{E l) (13) 
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depending on the Bose-Einstein distribution function 

ME) = ■ (14) 

It is well known, that this matrix constitutes a thermal Bogoliubov transformation, which 
combines particle and hole states in a way that guarantees a causal propagation in a thermal 
system 00. 

Historically, this transformation was discovered in the attempt to simplify perturbative 
calculations at finite temperature || . Hence, in the present section we are progressing in the 
reverse direction: We started from a simple perturbation expansion for the doubled fields, 
and then recover a causal description in terms of a Bogoliubov transformation of these fields. 

By comparison with standard results of thermo field dynamics, we obtain the Bogoliubov 



symmetry as |TC 



(E,k) (E,k)- jt ] (E, k) j W ( Ejk) = 

j{-\E, k) J{ + \E, k) - jt\E, k) 4 + \E, k) . (15) 

Of course, this relation admits many solutions for the Bogoliubov transformation matrix; 
it may be any element from the symplectic group in two dimensions Sp(2). Even the re- 
quirement of the KMS relation ( |10D leaves us with some freedom for the transformation 
matrix. However, it may be shown that for equilibrium states this freedom in the choice of 
the Bogoliubov matrix is irrelevant since it is a kind of gauge fixing |L0] - whereas for non- 



equilibrium states (where the parameter n is no longer an equilibrium distribution function) 
the choice of eq. ([T^) is the simplest possibility. 

The final step in our treatment of the generating functional for the free thermal Green 
functions is to obtain it in closed form for the causal sources j\ and jV Obviously the result 
is 



Z \j1J2] = exp 



-~ jd A x d 4 y 3a, x D% b (x - y)j b>l 



(16) 



The matrix- valued propagator in the exponent is conveniently expressed in momentum space, 
it has matrix elements 



Dl\k) = 


A F (k) - 


2nin B (uj k ) 5(kl -u\) 


Dl\k) = 


-2iri5{kl 


- ul) sign(fco) n B (k ) 


Df{k) = 


-2ni 5(kl 


- ul) sign(/c ) (l + n B (fco)) 


Df{k) = 


-am*) - 


2mn B (uj k ) 5{kl - u 2 k ) , 



(17) 

where is the Feynman propagator 

A F (k) = {k 2 -u 2 k + ie)- 1 (18) 

and A F its complex conjugate. The matrix elements are the four well known thermal Green 
functions: D u is the causal propagator, D 12 = D K and D 21 = D > are the thermal Wightman 
functions and D 22 is the anti-causal propagator. 



5 



III. INTERACTING FIELDS IN EQUILIBRIUM STATES 



In the previous section we have shown, how a technically simple description in terms of 
retarded and advanced propagators may be translated into a causal description through the 
use of a thermal Bogoliubov transformation. In the present section we extend this formalism 
to systems of interacting particles. For simplicity we assume, that other fields have been 
integrated out, i.e., the interaction appears in the form of a non-local self interaction. In 
general, also polynomial self-interactions like e.g. a 4 -theory fall under our description. 

The interaction is given by an arbitrary interaction Hamiltonian, which we express in 
terms of a field operator in the interaction picture labeled <p x . In this Hamiltonian, the 
field operators are replaced by field amplitudes belonging to the two representations of the 
doubled Hilbert space: 



(19) 



The primes in this equation indicate the fact, that so far we have not stated the temporal 
boundary conditions for the fields entering here. These follow from the almost trivial fact 
that also the interacting Green functions have to obey the KMS condition, i.e., the primed 
fields must be replaced by differential operators, according to 



-id 



Zo[h,32\ 



(20) 



Consequently, the generating functional for Green functions in a thermal equilibrium state 
and containing the proper boundary conditions is given by 



Z \j1J2] = exp -i / d 4 xH i 



-18 -iS 



(21) 



where we have to use the final expression for Zq of eq. flTE|). This generating functional 
is well known, and has been used very often in the existing literature. It gives the correct 
Schwinger-Keldysh structure of the perturbation series. 

However, the virtue of the path integral representation is that it goes beyond perturbation 
theory. Thus, we think it useful to apply the diagonalization concept of the previous section 
also to the interacting case presented here. Naturally this task is much more tedious than 
for the case of free fields. We therefore restrict the discussion to the level which is most 
important for any application, i.e., to the level of the full two-point Green function of the 
system. 

To this end we first define the generating functional for irreducible vertex functions in 
the usual way, cf. ref. [[U], p. 476]: 



r [0i, 2 ] = log [Z [ji, j 2 ]] - i / d 4 x ( ji{x)<j)i(x) - j 2 (x)(f)2{x) 



(22) 



T is the effective action functional for the interacting fields; for these we introduce a doublet 
to simplify the notation: 



^ 1 



01 (x) 

02 (x) 



(23) 
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A power series expansion for the effective action functional then yields 



T = l -d t $ T x r 3 d t $ x - U> T X r 3 (- A + m 2 ) 

jd i yU(x-y)% + 0(^). (24) 

Here, II (x — y) is the 2x2 matrix valued one-particle irreducible self energy function of the 
system, and t 3 = diag(l, —1). In principle one may also obtain the higher order terms of 
such an expansion, but for the purpose of the present paper we restrict ourselves to the first 
nontrivial piece involving the self energy function II. 

Since we deal with an equilibrium system, we may assume the self energy function to 
depend only on the difference of the space-time arguments. Consequently, one may perform 
the Fourier transform into 

F(E, k) = U(E, kf (r 3 (E 2 - ul) - U(E, k)) $(E, k) . (25) 

In accordance with the previous sections, the doublet is related to another field H by a 
Bogoliubov transformation |9|JT0[ 

$(E,k) =(B(n B (E)))- l ~(E,k) 

*(E, kf r 3 = E(E, k) T r 3 B{n B {E)) . (26) 

The fields H are therefore those which diagonalize the effective action functional up to second 
order. 

f(E,k)= (27) 

l^p.vr f E 2 -u 2 k -U R (E,k) \ i-^/j-, ,\ 

-^(E,k)r 3 ^ &-<4-UA{E,k)) 

T1 R and T1 A are the retarded and advanced self energy functions in momentum space, they 
are mutually complex conjugate and analytical functions of the energy parameter in at least 
one half plane: 

Tl R ' A (E, k) = ReU(E, k) t ma(E, k) 

= fde g(g ' fc) . (28) 

J E-E'±ie K J 

They are also related to the components of the matrix-valued self energy function by the 
simple relations 

U R {E, k) = U U (E, k) + U 12 {E, k) 

Tl A (E,k) = Tl u (E,k) + Tl 21 (E,k) . (29) 

More details of the diagonalization procedure, which (at least to the level of two-point 
functions, as considered here) may be formulated also for non-equilibrium states, can be 
found in refs. ||9"1,|T0 . 
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We proceed by reversing the Legendre transform, i.e., we want to establish how the 
interaction affects the diagonalization of the generating functional for Green functions. This 
reversal is straightforward, and one obtains for the generating functional in terms of the 
decoupled currents J\ and J 2 

Z[Ji,J 2 ] = 
exp 



- 1 - Jd*xd 4 y (j l)X D R (x - y)J ltV - J 2 , x D A (x - y)J 2 , w ) 
+ 0(J 3 ) 



(30) 



The retarded and advanced propagator occurring here are, of course, not the free propaga- 
tors. Instead they are obtained from a spectral function A(E, k) 



D R ' A (E, k) = ReD(E, k) =f mA(E, k) 



dE 



, A(E',k) 



E-E'±ie 

which in turn is related to the self energy function components by 

a(E,k) 



(31) 



A(E,k) 



(E 2 -u 2 h - ReIT(£, k)) 2 + n 2 a 2 (E, k) 



(32) 



The final step of the diagonalization procedure for the interacting case then is to perform the 
Bogoliubov transformation to the physical source currents conforming to the KMS condition. 
Quite similar to eq. (|l^) we obtain 



Z [jnh] = exp 



-- ld 4 xd 4 yj a , x D a \x~y)j b , y + 0(j 



(33) 



but the matrix valued propagator here is (in momentum space) 



D (ab \E,k)= JdE'A{E',k) 



x 



(34) 



(B(n B (E'))) 



( I 

V 



E - E 1 + ie 



B(n B (E'))r 3 . 



E-E'-ie/ 



Consequently, the thermal Bogoliubov transformation ( fl2|) diagonalizes the effective action 
functional up to the level of the full two-point function also for the interacting case. One may 
easily check, that in the limit of a free spectral function the full matrix valued propagator 
([M|) reduces to the free one, eq. (|17D . For the generating functional of Green functions, this 
diagonalization implies a factorization into upper and lower contour, see fig. 1. 
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IV. GENERATING FUNCTIONAL WITH INITIAL CORRELATIONS 



The next question to be addressed is how the above formulation has to be modified 
in the non-equilibrium case, i.e., how to treat particle propagation on a non-equilibrium 
background. To answer it we recall the principal feature of a non-equilibrium system: The 
occupation numbers of single particle states are not fixed a priory. For a perturbative 
treatment this implies some difficulty, since one may not be able to find operators which 
"annihilate" the non-equilibrium state. 

Consequently one may not use the "standard" Wick theorem to obtain the perturbative 
propagators. Necessary is a more advanced formulation of this important tool of perturba- 
tion theory. It has been shown some time ago that such a re-formulation amounts to the 
inclusion of initial correlations [ 17 , [I5| , , [I^| , which basically arise when a non-equilibrium 



"density matrix" is not a Gaussian in the full fields. 

The initial correlations appear in a perturbative expansion as additional vertex functions, 
coupling one, three, four, five, six or more fields, but not two fields. To explain this in 
more detail, we leave the path integral formulation for the moment, and resort to a Fock 
space picture with free field operators. To distinguish them from the C-number valued 
field amplitudes of the path integral formulation, we denote Fock space operators by a hat 
symbol, i.e., <f> x . The expectation value (•) denotes the trace over this Fock space with some 
non-equilibrium density matrix. 

In this picture we then apply Wicks theorem to a time-ordered product of L field opera- 
tors; resulting in a sum of products of normal ordered pieces and C-number valued contrac- 
tions. By normal ordering we will henceforth understand to shift all free-field annihilation 
operators to the rightmost side, i.e., a normal ordering with respect to the vacuum. However, 
our results also hold for more general normal ordering prescriptions with respect to an ar- 
bitrary equilibrium state: The argumentation is based on the undebatable fact that normal 
ordering is not possible with respect to a state that has unknown occupation numbers. 

Consequently, in a non-equilibrium state the expectation value of a time-ordered product 
of L free field operators is 

perm 71=0 

[L/2] 



n=0 



L ' CWiiD*)*, (35) 



(L - 2n)\n\2 



where ^ erm indicates the sum over all possible permutations of the L space-time coordinates, 
giving rise to the numerical factor in case of symmetrized space-time arguments (second line 
ofeq. (H)). 

Dll is the free causal propagator for the non-equilibrium state and the functions C^' 
are called correlation kernels . They may be expressed as a sum over expectation values 



of "vacuum-normal ordered" products. 

The advantage of the correlation kernels is that each of them appears only once in a 
specific order of the perturbative expansion. Note, that the only correlation kernel not 
present in a general non-equilibrium state is the two-point correlation: It is absorbed in the 
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propagator of the perturbative expansion. The odd correlation kernels only appear when 
the expectation values of the field unders consideration are nonzero. In this respect, the 
inclusion of initial correlations is a convenient way to obtain a perturbative expansion also 
for such fields without previous subtraction of the expectation value. 



DH is the 11-matrix element of 



D 



(oft) 
led 1 



x,y) = -i 




(36) 



A[-] denotes the anti time-ordered product. In equilibrium states this reduces to the matrix 
ofeq. (0). 

It is obvious, that the other matrix elements of this propagator occur in the Wick expan- 



sion of totally anti-time ordered, or mixed products (see ref. [ 1 1] for a complete discussion). 
It is therefore important to realize that the initial correlations do not distinguish between the 
two types of indices, i.e., they are not contour-ordered functions in the Schwinger-Keldysh 
formalism. 

We now drop the advantage of having a simple perturbative expansion, and replace 
the correlation kernels by their expansion into expectation values of simple field operator 
products. It has been shown in ref. | JL1|| , that this is also possible in terms of the Wightman 
functions rather than in terms of normal ordered products. This has the advantage that our 
result becomes independent of the type of normal ordering prescription we chose. Indeed, 
using for the unordered (Wightman) functions the abbreviation 



p{L) 
I...L 



(37) 



we may write 



[L/2] 



L! 



n=0 



(L — 2n)\ n\ (-2)' 



p(L-2n) ( / ( ' 



(38) 



Again, the numerical factor is due to the sum over all permutations of the space-time argu- 
ments, i.e., we understand the above expression to be symmetrized in these. 

The reason for this rewriting is, that one may easily obtain the generating functional for 
the correlation kernels in terms of the initial state Wightman functions, 



W [q] = exp ( +^ / d 4 xd A y i] x 



/Vy ) \^exp ( i j (I z //. ; 



(39) 



The correlation kernels are obtained from this by subsequent derivatives, i.e. 

'-i) K 5 K 



CI 



(K) 



,K 



W[rj\ 



The generating functional of the time-ordered functions therefore is 



(40) 
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exp ( i / d 4 z rj 



exp ( -- jd xd yr] x (iD nc0 (x,y) - {<p x <p v ) j //,, 



x ( exp ( i / d 4 z r] 



exp 



1 



d xd yr] x iD nc0 (x, y) x] y 



x W[rj] 



(41) 



From these expressions it is obvious how the equilibrium nature of a state is connected to 
the vanishing of the correlation kernels C: It is a gaussian of the fields 0, and therefore the 
generating functional W is identical to 1. 

Finally, we extend this formula to the matrix valued formulation. As we have indicated 
above, this brings into play the other matrix elements of eq. (Bfj). Furthermore the correla- 
tion kernels are identical for both types of indices, thus one has to replace the current rj in 
the functional W[rf\ by the sum of j\ and j%. The interaction is introduced in a perturbative 
way as in eq. (pl|) . The generating functional for non-equilibrium Green functions therefore 
is 



Zne [ji, j 2 ] = exp -i / d^u H { 



-iS -iS 



X 



exp 
exp 



- jd A xd A yj a ,x (iDf c0 {x,y) 



i / d 4 z (j hz + j 2 



X 



(42) 



This allows to recover the full perturbation expansion for non-equilibrium Green functions 
in an arbitrary state. States with nonzero field expectation value as well as states with 
complicated correlations are handled correctly by this expression. 

However, missing in this framework is the connection to the preceding sections, i.e., to 
the question of temporal boundary conditions. This gap will be closed in the following 
section. 



V. CUMULANT EXPANSION AND PATH INTEGRALS 



The expression for generating functional of free non-equilibrium Green functions in an 
arbitrary system is given in Eq. (^). We now proceed along the line of thought underlying 
the path integral formulation: Fields are expanded around their classical value, i.e., we set 
4> x = (fi x + 4> x . The classical field cl is the solution of the field equation with an external 
current as source 

(□a. + m 2 ) 4>t — J x ■ (43) 
Inserting this equation into (0) allows to rewrite the generating functional as 
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z [J] 



v\ 



exp 



- id 4 x(r x ( a x + ™ 2 )r x ) 



x exp 



- ld A x((f)p x ) 



(44) 



The first factor is a trivial Gaussian integral, i.e., it amounts to a phase of the whole 
expression and will be set to 1 henceforth. For cl in the second factor the classical solution 
obeying the correct boundary conditions must be substituted. In the case of equilibrium 
systems we may split up the fields and currents due to their location on the contour C 
depicted in fig. 1: 



d A y 



D R (x — y)J\, y if x on upper contour 
D {x — y)J2,y if % on lower contour. 



(45) 



This reproduces the generating functional from eq. fl5|). Another possibility for this splitting 
is the causal boundary condition 



d A y 



D (x — y)j a ,y if x on upper contour 
D 2a (x — y)j a ,y if x on lower contour, 



(46) 



which results in the expression ( [IB] ) for the generating functional. 

Here therefore is the place, where the temporal boundary condition enters the calculation 
of Z[j]. We now ask the question, how this classical equation has to be modified for non- 
equilibrium states. 

To this end, we introduce the irreducible cumulants by relating them to the generating 
functional of the correlation kernels ||19|| : 



d 4 x (j^Vx + Yl d 4 l...d 4 n 1 -C^ 

n=3 J U ' 



r)x . . . r} n 



exp 



d 4 x (4>Jj + d 4 l . . .d 4 n^ k(1, ...n)^ 

n=3 



Vn 



(47) 



In other words, the logarithm of the generating functional for cumulants equals the gener- 
ating functional for the correlation kernels. 

These cumulants n(l,...n) appear more than once in a given order of perturbation 
theory, whereas the correlation kernels appear only once. The use of cumulants therefore 
allows a resummation of the perturbation series similar to the Dyson equation; as well as 



the simple formulation of a linked cluster theorem |05 



For completeness we quote the expectation value of a time-ordered product in terms of 
these cumulants, i.e., the equivalent of eq. (p5|) : 



[L/2] L-2n 

EE E ^ 

perm n=0 j=0 



L-2n^ 



(48) 
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where K{(j) L }j denotes a product of j cumulants not containing 2-point cumulants whose 
arguments are a partition of the L arguments into j disjoint sets. 

Inserting these definitions in eq. (f|l|), we obtain for the causal generating functional in 
terms of the cumulants 



Z ne0 [ji, 32} = exp 



-~ JdVy f x iD: b c0 (x,y)j b y 
+ i Jd 4 x (^Jj (ji, x +j 2 , x ) 

00 „ . n n 

+ J2 ldH...d 4 n 1 -K(l,...n) HUx,i+ji, 



n=3 



(49) 



where = ji(xj), 2 = 1 . . . n. Note, that one may subtract the cumulant k(x, y) — \ ^x^y 
from the propagator and then obtains a complete series over all values of n for the last term. 
Also, the introduction of an interaction may be carried out as presented in section 3. 

We now use the results of section 2 to make the transformation into the re- 
tarded/advanced picture involving the currents J x . Although in section 2 we deal with 
an equilibrium state, it is clear that there exists a similar diagonalization transformation 
for the non-equilibrium case. To see this more clearly, we perform the Fourier transform of 
eq. (|36| ) with respect to (x — y) — > P, and write it in terms of the two arguments P and 
X = (x + y)/2. It follows that 



B{N XP ) (£>(«*)) r 3 {B{N XP ))- 1 



D* P (N XP D% P -(l + N X p)D% P ) 
D XP 



(50) 



where we may chose N XP in such a way that also the upper right corner of the propagator 
matrix is zero. In ref. |10| is is shown, that the corresponding N XP fulfills a transport 
equation and in the equilibrium limit approaches a Bose-Einstein distribution function. It 
is therefore indeed the generalized non-equilibrium distribution function. 

However, at the level of the external sources j, J this transformation has the disadvantage 
to be more complicated than eq. (|T2|). It mixes various momentum components of the ji to 
obtain the and vice versa ||10|| : 



d 4 xd 4 y 3a , x D( ab \x,y)j b , y 



d A q d A p 



(51) 



(2. )M 2./ 4Xe ~ l(9 ~ P)X (Si ) -3(TO(p + g) / 2 )))-V 3 



D (X, (p + ?)/2) 











D A (X, (p + q)/2) 



B(N(X, (p + q)/2))(f 



J2,q 



The "new" currents are marked by the horizontal braces. Let us note at this point, that the 
momentum mixing Bogoliubov transformation naturally couples the system to gradients in 
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the generalized distribution function. This allows to calculate transport coefficients [|1(J as 
well as non-equilibrium effects in two-particle interferometry |20 . 



This separation however fails when taking into account the higher order terms, i.e., those 
with correlations built in. Consequently, for a general non-equilibrium state the generating 
functional for Green functions does not factorize into a part of the upper branch and one 
on the lower branch of fig.l. This may also be seen from the equation that is fulfilled by the 
classical field 6 cl : 



Wit 

oo 

y- 



n-1 



d 4 l . . . d 4 n J\ . . . J n 



-(n-1)! 

1 c 

- Jd 4 yD la (x,y) Ja , y 

00 .„ „ n 

+ d 4 l...d 4 n Y[{ji,i + j2,i)«(l, ...,n,x). 



(52) 



n=3 



This equation contains the boundary condition for the classical field cl . Its is also ob- 
tained by application of a linked cluster theorem to the perturbation expansion, i.e., when 
evaluating the perturbation series for the classical field amplitude by hand considering only 
contributions from connected diagrams. 

Furthermore, it is obvious that it also approaches the Kubo-Martin-Schwinger (KMS) 
boundary condition (periodicity in the imaginary time direction) in the limit of zero cumu- 
lants, because it is then equivalent to eq. fl46|). 



VI. CONCLUSION 

With the present paper we achieved two goals. First of all we have extended the concept 
of diagonalized matrix valued propagators into the path integral formulation of quantum field 
theory. We were able to show, how to construct a simple perturbative expansion involving 
only retarded and advanced Green functions. This diagonalization procedure exploits a 
thermal Bogoliubov symmetry of the system || . The presence of this symmetry allows even 
simpler thermal Feynman rules when switching to a lower dimensional representation of the 
symplectic group fllPf . 

The diagonalization of the action integral amounts to a factorization of the generating 
functional for Green functions. Our result therefore includes the one of ref. |7[], but is more 
general since we have achieved the factorization also for the full two-point function. 

The second major goal we have achieved with the present paper is the inclusion of initial 
correlations into the path integral framework. This might be a future asset for two very 
different reasons. On one hand it allows to obtain a perturbative expansion also for truly 
non-equilibrium systems, where the boundary conditions for particle propagation are fuzzy 
at best. To this end one may use the results of section 4, where the correlations were 
introduced via the initial state Wightman functions. In particular, one may think of a 
practical application to transport theory, where the perturbative treatment only works for 
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some time. Thereafter one has to requantize the system with new degrees of freedom, e.g., 
with new quasi-particles. For the next time step then the correlations of the initial state 



play the role of additional vertices [IT . 

On the other hand, the correlations were related to a cumulant expansion of the action 
integral [TjJ . In this formulation, as presented in the last section of this paper, we were able 
to write down the integral equation fulfilled by the classical field amplitude. It is this integral 
equation which proves the close relationship between initial correlations and the boundary 
condition in time that one has to implement in quantum field theory for non-equilibrium 
states. 

We have furthermore proven, that the path integral for quantum fields in non-equilibrium 
states does not factorize into pieces on the upper and lower Keldysh contour. 
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FIG. 1. Time path for the action integral entering the generating functional. 
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